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We show how the rotational quantum state of a linear or symmetric top rotor can be recon- 
structed from finite time observations of the polar angular distribution under certain conditions. 
The presented tomographic method can reconstruct the complete rotational quantum state in many 
non-adiabatic alignment experiments. Our analysis applies for measurement data available with 
existing measurement techniques. 
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PACS numbers: 03.65.Wj, 33.15.Mt 

I. INTRODUCTION 

Recently, there has been great interest in non-adiabatic 
alignment of molècules using short non-resonant làser 
pulses , for a recent review see . Here the molècules 
are excited in a rotational state showing time-dependent 
angular anisotropy with ensuing revival structure long af- 
ter the passage of the aligning pulsc. The time-dependent 
angular distribution after passage of the làser pulse can 
be measured 0, 0, and the states created have many 
uses in ultra- fast òptics, high harmònic generation, scat- 
tering theory and potentially in chemical investigations 

So far, the comparison of such experiments with theory 
has been though the resemblance of experimental data 
with numerical simulations of the time evolution from 
a known initial state. Conversely, we propose a tomo- 
graphic method by which one can find the initial quan- 
tum rotational state of a linear rotor from measurements 
of the angular distributions at several points of time. 

A typical example of an alignment experiment is shown 
in figure ^ Here, a rotationally cold supersònic molec- 
ular beam is used as the molecular system in a pump- 
probe experiment. The initial state is thus approximately 
a thermal state with only a few energy levels excited. 
The pump is a short, with a typical duration of a few 
picoseconds, intense non-resonant làser pulse, which in- 
duces the rotational state. The pump pulse is linearly 
polarized and since the initial state is thermal and there- 
fore rotationally symmetric, the rotational state formed 
will be cylindrically symmetric around the pump polar- 
ization axis. Furthcrmorc, the field mainly interacts with 
the molecule through the molecular polarizability, mak- 
ing the rotational state reflection symmetric in a plane 
orthogonal to the pump polarization. This last symme- 
try is not, however, a requirement for the tomographic 
method below. 

After a delay, the molècules are coulomb-exploded by 
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Figure 1: The figure shows a typical alignment experiment 
with short làser pulses. The cold thermal molècules are ex- 
cited in a rotational state by a non-resonant pump pulse ver- 
tically polarized. After a time t the molècules are coulomb- 
exploded with an intense, ultrashort probe pulse, also verti- 
cally polarized. The fragments recoil back-to-back and are 
imaged on a 2-D detector. The 3-D angular distributions can 
be found from these 2-D images because of the cylindrical 
symmetry. Four examples of experimental data with Iodoben- 
zene are shown, where the imaged fragments are iodine ions: 
(from left) Probe applied long before pump, probe applied 
after respectively 1.5, 4 and 9ps. Part of the anisotropy, 
which is directly apparent in the un-aligned case, is due to 
the coulomb-explosion preferentially oceurring paral·lel to the 
probe polarization. This is corrected for when finding the 3-D 
distribution. Strictly speaking, Iodobenzene is an asymmetric 
top, but only slightly. For many purposes, it behaves like a 
symmetric top. Experimental data kindly provided by Simón 
S. Viftrup. 



an intense ultrashort probe-pulse, polarized paral·lel to 
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the pump pulse, and with a duration of typically some 
tens of femtoseconds. The molecular fragments recoil ap- 
proximately back-to-back and their angular distribution 
can be found from the imaged fragments. Because of the 
cylindrical symmetry, the true angular distribution can 
be found even though the fragments are pulled by an elèc- 
tric field towards a detection plane containing the làser 
polarizations. It is the purpose of this paper to present a 
tomographic method to reconstruct the rotational quan- 
tum state of the system from the measurement data in 
such experiments. 

We arrange the paper as follows: In section [l·l]we give 
a brief introduction to quantum state tomography. In 
section IIIII we present the tomographic method for the 
rigid linear rotor. In section ITVI we show how to extend 
the tomographic method to molècules where centrifugal 
distortion must be taken into account. In section Ivl we 
extend the tomographic method to the case of symmetric 
top molècules. In section fV*fl we discuss the possibilities 
of recording the experimental data necessary to perform 
a reconstruction and we conclude the paper. 



II. QUANTUM STATE TOMOGRAPHY 

Mathematically, tomography is a technique, closely re- 
lated to Fourier transformation, by which one can find 
a function from knowing its projection along all rotated 
lines. It found its first application in physiology in the 
1970s, and tomographic techniques are now used exten- 
sively in hospital scanners, where images of internal tis- 
sue can be found from data recorded using X-rays or 
NMR. Apart from physiology tomography is currently 
also being used in many other fields of research, includ- 
ing quantum physics and chemistry. Here the aim is to 
use measurements to find the quantum state of a system. 
The measurements are often taken to be spatial proba- 
bility distributions at different times, which for harmònic 
oscillator systems makes the method very similar to the 
scanning methods from physiology |(J, 0. Quantum 
state reconstruction has also been considered for many 
other systems including partides in traps (neutral àtoms 
0] and ions Q), general one-dimensional systems [ïoj| . 
dissociating molècules ^ÏJ and the angular state of an 
electron in a Hydrogen atom with n = 3 12]. 



III. RECONSTRUCTION METHOD FOR 
LINEAR MOLÈCULES 

In this section we will present a reconstruction method 
by which one can find the rotational state of a linear rigid 
rotor. As usually in quantum tomography we imagino 
that we are able to perform measurements of spatial dis- 
tributions at different points of time and that we know 
the Hamiltonian governing the time evolution. The goal 
shall be to find the unique quantum state corresponding 
to these measurements. 



A. Free rotation 

Before discussing the description of the quantum state, 
we shall first digress to account for what we know about 
the Hamiltonian and its eigenstates. As discussed in the 
introduction, we are varying the time elapsed between 
the rotational state is created by a pump pulse and it is 
probed. During this time the molecule is in a field-free 
environment. Consequently, we shall be interested in the 
Hamiltonian for the free rotor. 

The Hamiltonian governing the time evolution for the 
linear rotor with moment of inertia X is: 



H 



2J' 



where J is the angular momentum operator. The energy 
eigenstates of this Hamiltonian are { | J, m) } where the an- 
gular momentum quantum number J 6 No and the pro- 
jection quantum number m G {— J, — J+l, . . . , J}. These 
states have energy Ej — KIJ{J +1), where íl = £=. 
Since we will attempt to reconstruct the quantum state 
using angular distribution measurements only, we shall 
need the spatial coordinate representation of the energy 
eigenstates. As discussed further below, we shall restrict 
our study to azimuthally symmetric distributions, and 
we shall hence be interested in the distribution of the 
polar angle 9. For notational simplicity, we shall use the 
parameter x = cos(#) so the position representation of 
the eigenstates becomes: 



t(x\J,m) = P?(x)e 



-iQJ(J+l)t 



(D 



There should be a factor of e %ra ^ in Eq. JIJ, but we can 
safely ignore this since we shall restrict our treatment to 
fixed vàlues of m [24j • The PJ 1 are the normalized asso- 
ciated Legendre polynomials, normalized on i £ [-1, 1] 
HI: 



dxPf x {x)P^{x) = 6 JuJa . 



(2) 



The multiplicative factors relating these to the un- 
normalized associated Legendre polynomials PJ L {x) are: 



'2J+1 (J-m)! 



(J + m)! 



P?(x). 



Eq. 0} fully accounts for the free time evolution of the 
rotor. The angular position distribution of a single eigen- 
state is constant in time, but due to the different time 
dependent phase factors, a linear combination of differ- 
ent | J, m) states will have time dependent interference 
terms. Hence, such a linear combination will have a time 
dependent position distribution, also in the field free case. 

After this brief account of the dynamics of the system, 
we proceed to define what we mean by the quantum state 
we are trying to reconstruct. 
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In the present treatment, the state will be character- 
ized by its density operator p, which allows a descrip- 
tion of statistically mixed states. Specifically, we shall 
find all the matrix elements of p in the | J, m)-basis at 
t = restricted to single m's. This is the same as finding 
the diagonal blocks of constant m in the density matrix, 
i.e. p(J,m, J' , m) = p m (J,J'). For a general state this 
is only a partial characterization, but if it is somehow 
known that the state has no correlations between differ- 
ent m vàlues, knowledge of all p m {J,J') amounts to a 
full characterization of the state. This is the case in the 
experiment described in the introduction and the proce- 
dure below will completely reconstruct the quantum state 
for such systems. For the angular distribution, this di- 
agonality in m implies cylindrical symmetry for all time 
around a laboratory fixed axis, which we choose to be a 
spherical polar axis. 



B. Measurements and observables 

Having described the physical system, we proceed to 
consider the measurements performed; namely the (an- 
gular) position distribution measurements at the time t: 



Writing the position distributions Eq. J3J in the new vari- 
ables J and AJ, and using Eq. we arrive at: 



Pr(x,t) 



oo J oo J 

Y Y + Y Y 

V 



J=\m\ AJ=-J J=\m\ AJ=-J . 
J even AJ even j dd AJ odd/ 



fJ + AJ J-AJ 
e -inAJ( J+ x)t J2 Vl{x)Cr AJ , L . (5) 

L=\AJ\ 

For readability, we shall in the following abbreviate the 
summations in the first line of Eq. (J5J as: 



E- 

J, AJ 



(6) 



We have now expanded the time dependent position dis- 
tribution on orthogonal polynomials V^(x) with coef- 
ficieents that can be found by simple spatial integrals 
J^ 1 dxV™(x)Pr(x,t). These coefficients, in turn, are lin- 
ear combinations of the density matrix elements which 
can be found if the distribution is recorded at different 
times, as we show below. 



Pr(ar,í) = t (x\p m \x) t 



OO OO 



= ^ ^ t{x\Ji,m){Jum\p m \J2,m){J2,rn\x)t 

Ji — \m\ J 2 — |m| 



oo oo 



= E E PMi,J2)VT 1 (x)VT 2 (x)x 

Ji — |m| J 2 — |m| 

e -ií2[./i(Ji + l)-./ 2 (./2 + l)]t (3) 

where t = is an arbitrary, but fixed, point of time af- 
ter passage of the pump pulse. Our aim is to invert this 
equation to find the density matrix in the eigenstate- 
representation p m (Ji, Ji). It will prové convenient to 
write the produets of the two Legendre polynomials 
Vj l (x)Vj^(x) as a sum of single Legendre polynomials. 
This is exactly what is done in the decomposition of di- 
rect product bases for irreducible representations of the 
rotation group: 



J1+J2 

= Ev 

L = \.h-J 2 



2%C(J 1 ,J 2 ,L\m,-m,0)Vl{x), (4) 



where C( Ji , J2, L\m, — m, 0) are Clebsch-Gordan coeffi- 
cients [l4|. Next, we introduce the variables J = J\ + J2 
and AJ = J\ — J2, and the notational simplification: 



C. Reconstruction formulas 

In this subsection we will present the reconstruction 
formulas. We will first consider the off-diagonal elements 
of the density matrix and deal with the diagonal clements 
further below. 



1. Reconstruction of off-diagonal elements 

We let NtT be the time interval in which position 
measurements Pi(x, t) have been performed. Here Nt € 
N and the minimum measurement time T = 2ir2/h, 
equalling \f2 times the semiclassical rotational period of 
the J = 1 state or the rotational revival period ■ 

The time dependent factor e ^ l ^ AJ ( J + 1 ) t m Eq. @ sug- 
gests that a temporal Fourier transform evaluated at fre- 
queney £lf3(a + 1), with a and (3 suitably chosen integers, 
will be part of the reconstruction procedure. Following 
this strategy, we consider the integral: 



1 



NtT 



1{a ^ ) = làr Jo dteÜma+1)t JdxV° a (x)Pr(x,t) 
ÍJ + AJ J-AJ\ 

2 , Pm\ 5 ' 9 / Ó / 3 ( Q + 1 )- AJ ( J + 1 ) X 

J,AJ ^ ' 



c 



J,AJ-L 



J + AJJ-AJ 



m, —to, 



E C ™AJ,L S L,a, (7) 
L = |A,7| 
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where a S {\m\, \m+l\, . . .} and f3 S {— a, —a+2, . . . , a}, 
which makes both a and (3 are either even or odd. For 
the time being, we shall choose (3 ^ and deal with the 
(3 = case (i.e. diagonal matrix elements) later. 

Notice that the delta- function in the sum over L en- 
sures that a, J and A J are all of the same parity. This 
is also the parity of /?, since we chose a and (3 to be of 
the same parity. 



J=|m| A.J=-J 
parity as a parity as a 



<5/3(q + 1)-A./(J+1)C™ A j 

L=|A,/| 

Rather than giving a single term in the general case, 
the sums are greatly simplified. From the L-sum we get 
I(a, (3) = unless: 



|AJ| < q < J. 



(9) 



Using this in conjunction with the other delta-function 
we find I(a, (3) = unless: 

|AJ| < \f3\ < a, (3 and A J of same sign. (10) 

Furthermore, this delta-function demands: 

Í3(a + 1) = AJ(J+1). (11) 

As it turns out, for many choices of a and (3 there is 
only the straightforward solution to the equation Eq. Ijlll) 
under the conditions Eqs. JÜJ) and l|10[) . namely a = J and 
j3 = A J. In these cases we easily find the density matrix 
elements: 



I(a,/3) = C£ ataPm 



^ [3 (X _ \ uniquc solution/ 1 rj\ 
2 ' 2 / ' tol9ltÏ0land[TïlV 1Z ) 



For some choices of a and (3 there will be more than 
one term surviving from the sums in Eq. (fHJ . The 
physical reason for this is that there is more than one 
pair of energy eigenstates having a certain energy dif- 
ference. One example is the energy difference between 
the states (Ji = 3 — > J2 — 0) and between the states 
(Ji = 6 — > J2 — 5). In our variables this corresponds to 
(J = 3, AJ = 3) and (J = 11, AJ = 1). We proceed to 
show that all density matrix elements with (3 ^ can be 
found regardless of this complication. 

Imagine that we have found I(a, (3) to contain 

the term C 



temis CE a , a 



( a+p a-/3 \ 
l \ 2 ' 2 I ' 

^m\^ 2 ' 2 y 



but also several other 



where 



"1 



1, 2, . . . , JVi. The strategy is to calculate all the integrals 
I(a' n ,/3' ). These may also contain several terms with 
(J = < 2 , AJ = ffj, where n 2 = 1, 2, . . . , N 2 , but al- 
ways fewer than in the former integrals, i.e. Nj < Nj + \. 
The reason for this is the condition |AJ| < \(3\ from 



Eq. (|10() : Whereas the (3 = AJ term will always appear, 
the rest of the terms in the sum Eq. © will have A J < /3, 
whereby the procedure will tcrminate and we can find all 
the matrix elements of p m by back substitution. If one 
is interested only in matrix elements up to some max- 
imum J = J ma x it is only necessary to search for de- 
generades in the energy differences for all (J, (3) up to 

J Jrnax{.Jmax f)/2 for even Jmax and Jraax{,Jmax+l) 
for odd J m ax- 

To clarify this procedure a little, let us consider 
an example: We shall choose m = and try to 
find the matrix element po(5,0). In this case J = 
AJ = 5, and we find 7(5,5) = C^ 5 5 p m (5, 0) + 
C 9?3, 5 P™(M) + C^ 15 p m (15,14). To find 'the desired 
matrix element we shall need the two extra equations 
7(9,3) = C^ 3 pm(6,3) + C^ h3 p m (15,U) and 7(29, 1) - 
C29 1 iPm(15, 14). These three resulting equations can 
now be solved by back substitution to yield /? TO (5,0), 
p m (6, 3) and p m (15, 14). 



2. Reconstructing the diagonal 

We proceed to show how one can go about finding the 
diagonal of p m in the energy-representation. To do this 
we shall use Eq. © with j3 = 0, by which a is even: 



l(a,0) 



1 



00 

J=\m\ 
J even 



N T T pi 

dt JdxV°(x)Pr(x,t) 



J J 

2' 2 



C 



J,0,a 



Pm{ J\, Jl) C™ Jl al 

Ji — \m\ 



(13) 



where we have used J\ = J/2. We notice that the coef- 
ficient C^} 1 a is zero unless a < 2Ji. We arrange the 
J(a, 0) in a column vector 7 and the p m (Ji, Ji) in a col- 
umn vector R: 



( /[2|m|,0] 
Z[2(H + 1),0] 



p m (|m|, |m|) 
_r = I Pm{\m\ + 1, jmj + 1) 



Eq. I|13|) now reads: 



<ét ■ R, where 

I °2|m|,0,2|m| L· '2(|m| + l),0,2|m| 



(14) 



a 



2(|m| + l),0,2(|m| + l) 



upper— triangular 
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To find the desired elements p m (Ji, Ji) we shall need to 
invert the matrix „#, which can be done since the diag- 



onal elements C™ a = y / 2nC(a/2, a/2, a\m, — m, 0) are 
all non-zero. In practical applications, one may truncate 
the linear system at some maximum J\ corresponding to 
a maximum energy. 

It may be noted that since the linear system (|14|) al- 
ways has the same coeficient matrix, we could just has 
well have used a function J-j 1 that selects a certain Jj in 
(|T3j) . so that: 

r N T T f l 

pM,Ji) NtT 

The Tj x (x) can be precalculated and are given by: 



— dt dxF^{x)??{x,t), 

T± Jo J-l 



OO 

E 

J=Jl 



fj u jV2j(x), 



where the expansion coeficients fj u j are easily found by 
Cramer's rule and the fact that ^# is upper-triangular: 



ÍJuJ = (-1) 



Ji+J 



'-Ji,J\ 



and 



ysub 
l Jl,J 



?yfj | is the determinant of the sub-matrix: 
Jl (Ji + 1 : Ji + 1 + J, Jx + 2 : Ji + 2 - 



(15) 



This mcthod would bc relevant if only a few particular 
elements of the density matrix are desired. One may 
notice the similarity of these functions to the pattern 
functions used with harmònic oscillator systems [l5 | .[Ï6 j . 

It may be noted that we could not have circumvented 
the back substitution procedure in section IHI C II by 
choosing other functions than the V® (x) in Eq. © , since 
the solution of the resulting system of equations depends 
on the measurement vàlues through the /(<x/, /?j)'s. 



IV. REAL MOLÈCULES AND CENTRÍFUG AL 
DISTORTION 

In the above treatment we used the rigid rotor Hamil- 



tonian H 



21 ' 



whereas terms of higher order in J will 



be important in real molècules in highly excited rota- 
tional states. With the lowest order centrifugal distortion 
the Hamiltonian is instead: 



H' 



21 
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The states \J, m) are still the energy eigenstates, but here 
with energy Ej = MÏJ(J + 1) - HDJ 2 (J + l) 2 . The 
constant D oc j$ ± , where u> is the harmònic frequency 
of the molecular bond. The reconstruction strategy is 
similar to the one in Eq. JJJ: 



lim 



T 

' dte if}(a+l)[n-Df3( a +l)]t x 


1 

dxVl(x)Pr(x,t). (16) 



In this case, the conditions Eqs. JÏÏJ and l|10|l still apply, 
but Eq. 111(1 is replaced by: 



AJ(J + 1)-^AJ 2 (J + l) 2 



(17) 



Depending on the value of D/Çl, there may be more than 
one set of (A J, J) that satisfies the conditions Eqs. ©- 
HI U| l and Eq. I|17|l . but it will usually be much fewer 
than in the rigid rotor case. Consequently, one can of- 
ten avoid the back substitution procedure above and use 



V. SYMMETRIC TOP STATES 

Having treated the linear rotor case, we will now 
straightforwardly generalize the results from section IIIII 
to symmetric top molècules. This class contains many 
more molècules than the linear class. 

The eigenvectors of the Hamiltonian are in this case 
the set {| Jkm)} with associated energy E = híl\J(J + 
\) — hQ,2k 2 . Here k is the quantum number for the projec- 
tion of the angular momcntum on the symmetry axis of 
the molecule, where m was associated with a projection 
on a space-fixed axis. Like m, k can assume the vàlues 
— J, — J + 1, . . . , J. For a linear molecule k = 0, and the 
treatment in section IIIII can indeed be seen as a special 
case of the symmetric top. In accordance with the above 
treatment we shall assume that we can measure the polar 
angular distribution, parametrized by x = cos(0), for a 
certain value of k and m. The angular position represen- 
tations of the eigenstates are: 



t(x\ Jkm) 



2J+1 



d J km {x)e 



-í(fii j(j+i)-n 2 fe 2 ) 



where the d/, m (x) are the usual rotation matrix elements, 
and the eigenstates are normalized on [—1, 1]. The quan- 
tities íli = l/2X x and ÍI2 — h/(2T z — 2I X ). Here we 
have chosen the z axis as the molecular symmetry axis, 
and used that the moments of incrtia T x = I y . As in 
Eq. we can safely ignore the dependence of the two 
other Euler angles <$> and Xj since we are working in the 
subspace of fixed k and m. We can now readily generalize 
the treatment in section IÏTÏ1 bv exchanging the functions 
V'j l (x) with y/ (2 J + l)/2 dj£ m (x) and the lower summa- 
tion indices \m\ with M krn = max(|m|, \ k\). In particular, 
Eq. generalizes to: 



Pr(x,í) = t(x\pk, m \x)t 

OO OO 

=E J2Pm(^^)dL^)d J k 2 m (x) 

Jl= J2= 



(2Ji + 1)(2J 2 + 1) iU[ j l( j 1+1 )_j 2( j 2+1 )] t 
4 



(18) 
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Products of the functions d^. m {x) can again be decom- 
posed into linear combinations of single Legendre poly- 
nomials, again introducing J — (Ji + J 2 )/2 and AJ = 
(Ji - J 2 )/2 0, 0: 



(2Ji + 1)(2J 2 + 1) 



t=|Ji-J 3 | 



(19) 



where 



2L + 1 



C( Ji, J2, L\m, —m, 0) x 
C(J u J 2 ,L;k,-k,0). (20) 



r~ikra 



Proceeding as in section IIIII we calculate the integrals 
I(a, (3) and arrive at the equivalent of Eq. © : 

/(a,/3) =2^ 2^ pk, m l , 1 x 

J=M km AJ=-J ^ ' 

parity as a parity as a 

J 

àp(a+l)-AJ(J+l)Cj%J,a X/ ^ L ' a ' ( 21 ) 
L = \AJ\ 

The reconstruction of the matrix elements Pk,m{Ji, J2) 
can therefore be done like in sections IIII C II and IIII C 21 
where one uses the C k /\j l defined in Eq. 112011 instead of 



VI. DISCUSSION 

As the reader may have noticed, the method above 
effectively deals with a semi-continuous one-dimensional 
reconstruction in the space of x and J, and not with the 
reconstruction of a general state. The reason for this lim- 
itation is that the inversion formulas we are using are in- 
tegral transformations that preserve dimcnsionality, like 
e.g. the Fourier transformation in t in Eq. ©. If the 
dimensionality of the measurements are not at least as 
large as the dimensionality of the quantum state, such 
an approach nearly always fails |l9j |. For example, in one 
dimension the density matrix is a two-dimensional ob- 
ject (x\pk,m\%'} which has the same dimensionality as the 
measurement-space Pr(x,t). In contrast, the full quan- 
tum state of e.g. the linear rotor is a four-dimensional 
object (x, 4>\p\x' , (/>'), while the angular measurements 
Pr(x, <j), t) would only be three-dimensional. If one de- 
sired to reconstruct a full general quantum state, one 
would either have to be able to perform more advanced 
measurements or to vary the Hamiltonian as suggested in 
[20| . Because of the difficulty of precisely knowing Hamil- 
tonians arising from e.g. additional làser pulses, we have 
instead treated the one-dimensional system using only 
routinely performed measurements. 



In spite of this limitation, the situation where the one- 
dimensional description is a complete characterization of 
the quantum system is experimentally common. Indeed, 
for the situation described in the introduction, there are 
no correlations between the different k and m quantum 
numbers. The question is of course whether one can per- 
form measurements of polar angular distributions for cer- 
tain k and m quantum numbers. 

A conceptually straightforward solution would be to 
use very low-temperature samples, where only k = m = 
is populated before the interaction with the pump pulse. 
Note that this would not imply that the created rota- 
tional state would be a pure state due to inhomogencitics 
in the pump pulse focus, e.g. the pump intensity gener- 
ally varies over the focus or shot-to-shot variation of làser 
intensity. Though conceptually simple, this situation is 
unpractical due to the smallness of the rotational energies 
demanding very low temperatures. 

A more practical approach would be to select a cer- 
tain value of k and m before interaction with the pump 
pulse. For polar symmetric top molècules this can be 
done for certain vàlues of k and m by using the lin- 
ear Stark effect through hexapole focusing techniques 
[2Ï|. p^ . The linear Stark shift for symmetric tops is 
AE^ 1 ' = —pekm/J(J + 1), where p is the dipole mo- 
ment and e is the magnitude of the elèctric field. Sincc 
one starts out with a thermal sample, there are no corre- 
lations in either J, k or m in the state before interactions 
with the pump. Therefore, one does not discard any in- 
formation on the state by picking out all sets of vàlues of 
J, k and to one by one. If the molecular beam was suf- 
ficiently cold, one could select components of the beam 
with certain k and m vàlues, as is required by the method 
above. One would then also select a certain value of J, 
which is not required, but on the other hand poses no 
problem. Furthcrmore, the apparent ambiguity that a 
change of sign of both k and m gives the same Stark 
shift is of no consequence, since the state created by 
the pump l·icLS idcnticcil ma-trix clements pk,m 

(J, J') and 

P-k,-m(J, J') due to the invarianec of the initial state 
and of the pump-pulse interaction Hamiltonians under 
the interchange (fc,m) — > (— k, —to). 

Most linear molècules and symmetric tops with k = 
do not exhibit any linear Stark shift [2Íj, but also here 
hexapole focusing techniques can be used to select spe- 
cific vàlues of J and m through the second order Stark 
shift. Here, both molècules with the quantum numbers 
(J, to) and (J, — m) are selected, but this is of no conse- 
quence since p m (J,J') = p- m (J,J') in accord with the 
previous paragraph. 

In conclusion, we have presented a method to recon- 
struct the blocks Pk,m( J, J') of the rotational density ma- 
trix for linear and symmetric top molècules from angular 
distributions. This amounts to a complete characteri- 
zation of the quantum state in a common experimental 
setup. Finally, we have suggested how the required mea- 
surements could be performed with existing techniques. 
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